Topological superconductor to Anderson localization transition in one-dimensional 

incommensurate lattices 



(N 

o 

(N 
< 

m 

(N 



Xiaoming Cai, 1,2 Li-Jun Lang, 1 Shu Chen, 1 ^ and Yupeng Wang 1 



jing National Laboratory for Condensed Matter Physics, 
Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
'' Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, 

(Dated: August 24, 2012) 



China 



We study the competition of disorder and superconductivity for a one-dimensional p-wave super- 
conductor in incommensurate potentials. With the increase in the strength of the incommensurate 
potential, the system undergoes a transition from a topological superconducting phase to a topo- 
logically trivial localized phase. The phase boundary is determined from various aspects and the 
topological superconducting phase is characterized by the presence of Majorana edge fermions in the 
system with open boundary conditions. We also calculate the topological Z2 invariant of the bulk 
system and find it can be used to distinguish the different topological phases even for a disordered 
system. 
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Introduction.- Topological superconductors (TSCs) 
have attracted intense recent studies, as they are promis- 
ing candidates for the practical realization of Majo- 
rana fermions [IHtJ- Among various proposals, the one- 
dimensional (ID) TSC in nanowires with strong spin- 
orbit interactions and proximity-induced superconduc- 
tivity 0, 0] provides experimental feasibility on the de- 
tection of Majorana fermions in hybrid superconductor- 
semiconductor wires [il-fioj]. which has stimulated great 
enthusiasm in exploring physical properties of topologi- 
cal superconductors. A key feature of a ID TSC is the 
emergence of Majorana edge states at ends of the super- 
conducting (SC) wire as a result of bulk-boundary corre- 
spondence. A prototype model unveiling the topological 
features of the ID TSC is given by the effective spinless 
p-wave SC model studied originally by Kitaev [2|. 

As the TSC is protected by the particle-hole symmetry, 
the topological phase is expected to be immune to per- 
turbations of weak disorder 11] . Nevertheless, a strong 



disorder may destroy the SC phase and induce a transi- 
tion to the Anderson insulator. Localization in ID SC 
system in the presence of disorder has been an active re- 
search field in the past decades 12h15|. The theoretical 
studies have unveiled that the particle-hole symmetry in 
the SC system plays an important role in the problem 
of the Anderson localization 12|. Due to the existence 



of a finite SC gap, the interplay of disorder and super- 
conductivity leads to a topological phase transition from 
topological SC phase to a topologically trivial localized 
phase when the strength of disorder increases over a crit- 
ical value. 

So far, most theoretical work for the Anderson local- 
ization in ID TSCs focuses on the random disorder 
17] , disorder produced by incommensurate potentials is 
concerned only very recently [3]. In the absence of 
superconductivity, the localization transition driven by 
the incommensurate potential occurs at a finite disor- 



der strength which can be exactly determined by a self- 
duality mapping [l9| , whereas an arbitrary weak random 
disorder induces the Anderson localization in one dimen- 
sion. The incommensurate potential can now be engi- 
neered with ultracold atoms loaded in ID bichromatic op- 
tical lattices [20I ] , opening the experimental way to study 
the localization properties of quasi-periodic systems. In 
this work, we shall study the interplay of the incom- 
mensurate potential and topologically protected super- 
conductivity in the ID p-wave SC model and determine 
the phase boundary of TSC to localization transition in a 
numerically exact way. The tunability of the incommen- 
surate potential [20j j provides a potential way to exper- 
imentally study the controllable disorder effect in TSCs 
realizable in cold atom systems plj . 

Model of p-wave superconductor with incommensurate 
potential.- The ID p-wave superconductor in the incom- 
mensurate lattices is described by the following Hamilto- 



H.c.) + Vifn 



(1) 



where hi — c\ci is the particle number operator and 
cj (cj) the creation (annihilation) operator of fermions. 
Here the nearest-neighbor hopping amplitude t and the 
p-wave pairing amplitude A are taken as real constants, 
whereas the incommensurate potential 



Vi = V cos(27ria) 



(2) 



varies at each lattice site with a being an irrational num- 
ber and V the strength of the incommensurate poten- 
tial. The model reduces to the Aubry- Andre model when 
A = [13, while the Hamiltonian describes the Kitaev's 
p-wave SC model for a = [2]. For A = 0, the sys- 
tem undergoes a derealization to localization transition 
at V = 2t. On the other hand, the uniform p-wave SC 
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system with V i = V undergoes a topological phase tran- 
sition at |V| = 2t with a topological nontrivial phase in 
the regime of |V| < 2t characterized by the presence of 
edge Majorana modes Q. In this work, we shall study 
the interplay of the SC pairing A and the incommensu- 
rate potential and then determine the phase diagram of 
the system. 

The Hamiltonian can be diagonalized by using the 
Bogoliubove-de Gennes (BDG) transformation defined 

as, 



i=l 
L 

Vn = ^2i u nA \+V n ,A], 
i=l 



(3) 



where L is the number of lattice sites and n = 1, • • • , L. 
Here u n ^ and v n ^ are chosen real. In terms of the 
operators rj n and 77^, the diagonalized Hamiltonian is 

written as H = J^n=i ^niVnVn ~ 5) witn Ki being 
the spectrum of the single quasi-particles. The spec- 
trum as well and v n i can be determined via 
the diagonalization condition, [rj n ,H] = A n r] n , which 
leads to the following coupled equations: A n u n ,i = 
-t(u n<i+ i + itn.i-i) + ViU n .i - A(v n . i+1 - v n<i -x) and 

A„U„,i = t(v n ,i+i + Vn,i-l) - ViV n<i + A(u nt i +1 - U n ,i-l) 

with i = 1, • • • , L. The above equations can be also writ- 
ten as a more familiar form, i.e, the BDG equations, as 



h A 
-A -h 



= A r 



(4) 



where 



-t(5j,i+i + 8j t i-i) + ViSji, 



A 



1 j 



-A(<y Ai 



Here = (u n<1 ,--- ,u n<L ) and v% = (u„,i,--- ,V n ,L)- 
Given (tt n> i, u„,i) with i = 1, • ■ • , L being the solution to 
Eqs. (@} with a positive eigenvalue A„ > 0, one can prove 
that (v n< i, u n> i) is also the solution to BDG equtions with 
the eigenvalue — A„, which implies ?7 n (A n ) = A„). 
The ground state of the system corresponds to the state 
with all negative quasi-particle energy levels filled. If the 
quasi-particle energies are arranged in ascending order, 
i.e., Aj < Aj+i, for A, > 0, the gap of the system is just 
given by A g — 2Ai . In the following calculation, we shall 
set t = 1 as the energy unit. 

Transition from SC phase to disorder phase.- For the 
incommensurate potential with an irrational a, the so- 
lution to BDG equations equals to solving an eigenvalue 
problem of a 2L x 2L matrix. Solving Eqs. (0}, we can 
get the whole spectrum of quasi-particles. In Fig.l, we 
show the spectra for the case of a = — l)/2 and 
A = 0.5 under periodic boundary conditions (PBC). It 




FIG. 1: Energy spectra of ID p-wave superconductors with 
a = (x/5 - l)/2, A = 0.5 and L = 500 under PBC. 
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FIG. 2: (a). The energy gap_A 9 versus V - 2 A, (b). the 
average correlation function Cl/2 versus V — 2 A, (c). the 
MIPR versus V - 2A for the system with a + (\/5 - l)/2 and 
L = 500. 



is shown that there exists a regime with obvious nonzero 
gaps when V is smaller than a critical value V c . When 

V exceeds the critical value, there is no an obvious gap 
separating the negative and positive parts of spectra. To 
sec it more clearly, we show the variation of A g versus 

V in the regime close to the transition point in Fig. 2a. 
As shown in the figure, the gap vanishes at about V c = 3 
and the system opens a very narrow gap in the regime of 

V > V c - We calculate the gap for systems with different 
A and find the similar behavior: the gap reaches a min- 
imum, which approaches zero in the limit of L — > 00, at 
the transition point about V c = 2 + 2 A and there exits a 
very narrow gap when V exceeds the transition point. 

Observing that the p-wave fermion model corresponds 
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to the transverse XY model with a randomly (irra- 
tionally) modulated transverse field [22|, [23jj: H = 

~ J2i[ J x a i a i+i + Jy<Ti<7i+i] + J2i h i a i> with tne identifi- 
cation of J x = (i + A)/2, J y = (t-A)/2 and h { = -Vi/2, 
we can identify the phase transition by calculating the 
correlation function Cjj = (crfuj). In the language of 
quantum spin model, the ferromagnetic phase is charac- 
terized by the long-range order of the correlation func- 
tion (erf crj) >co = A with A being a nonzero posi- 
tive number. In the original fermion representation, erf — 
(cj + Ci)cxp(— in CjCj) takes a nonlocal form includ- 
ing a string product of fermion operators, and the corre- 
lation function Cy = ((q+Cj) exp(— wrj^^ n;)(cj+Cj)). 
In the presence of the disordered potential, the correla- 
tion function dj will oscillate and we define the average 
correlation function C r — J^i ^i.i+r/L. Then for a large 
system under PBC, the value of Cl/2 can be used to 
distinguish the SC phase and the localized phase. The 
correlation function Cy can be calculated by the exact 
numerical method described in Ref.[23]. In Fig. 2b we 
show the relation between Cl/2 and V for systems with 
different A. Without the disordered potential, the corre- 
lation function is a positive number and increases as 
A increases for < A < 1, gets its largest value Cl/2 = 1 
at A = 1 , then decreases for A > 1 . As the strength of V 
increases, Cl/2 decreases monotonically and approaches 
zero when V-2A is about 2. When V > 2 + 2 A, the sys- 
tem loses the long range order of correlation function and 
the system is driven into the Anderson localized phase. 

When V > V c , we would expect that there exists 
a transition from the SC phase to the disorder phase. 
To characterize the localization transition, we define the 
quantity of the inverse participation ratio (IPR) which is 
defined as 



Pn 
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FIG. 3: The finite size analysis of the transition point, i.e., 
Vc(L) vs 1/L. Inset: Ai(L) = A g (L)/2 vs 1/L in the regime 
of V < Vc and V > V c . 
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FIG. 4: Energy spectra of ID p-wave superconductors with 
a = (VE- l)/2, A = 0.5 and L = 500 under OBC. 



where u n j and v n j are the solution to the BDG equa- 
tions and fulfil the normalization condition X)i( u ni + 
v ni) = 1- The definition in Eq.© can be viewed as an 
extension of IPR for the case with A = [24 , 25 1 . For 
an extended state, P n — > 1/L and the IPR tends to zero 
for large L, whereas the IPR tends to a finite number 
for a localized state. Therefore, IPR can be taken as 
a criterion to distinguish the extended states from the 
localized ones. Since the ground state is composed of 
states with all negative quasi-particle energy levels filled, 
we define the mean inverse participation ratio (MIPR) as 
MIPR = J2n=i-Pn/L to characterize the localization of 
the ground state. As shown in Fig. 2c, the MIPR increases 
monotonically with the increase of V. At V = 2 + 2A, 
the MIPR has a sudden increase which characterizes a lo- 
calization transition. As a comparison, we note that the 
localization transition does not occur for the commensu- 
rate potential system with a rational a 26], for which 



nvi n ), we calculate the 
„,„ x and L m in and find 



the wave functions of a periodic system take the Bloch's 
form and are extended for arbitrary V . 

To determine the transition point more precisely in the 
limit of L — > oo, we make finite size analysis by calculat- 
ing the transition points for systems with different sizes. 
As shown in Fig. 3, the value of transition points V C (L) 
for systems with A = 0.5 oscillates around 3.0. Defining 

Vavc = L=L min Vc\L)l {Lraax ~ L 

average of V C (L) for different L 
that Vavc is about 3.0040 ± 0.0005 being very close to 3. 
The change of the gap size at V — 2.5 and V — 3.5 is 
shown in the inset of Fig. 3, which indicates that the gap 
is finite in the regime of V < V c whereas the narrow gap 
in the regime oiV>V c approaches zero in the large L 
limit. We also check systems with different A and find 
similar behaviors, i.e., V C (L) — 2A oscillates with L and 
approaches to 2.0 in the large L limit. 

Topological features of the topological SC phase.- To 
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FIG. 5: The spatial distributions of <f>i (a) and ipi (b) for the 
lowest excitation with various V . 



characterize the topological properties of the SC phase, 
one need seek the zero-mode solution of the system un- 
der open boundary conditions (OBC). As shown in Fig. 4, 
we plot the quasi-particle spectra of BDG equations un- 
der OBC. In comparison with the spectra under PBC, 
an obvious feature is the the presence of the zero mode 
solution in the gap regime. To see clearly how the zero 
mode solution evolves with the change of the disorder 
strength, we show the enlarged Ar in the inset of Fig. 4. 
It is shown that there is a sudden increase in Ai when 
V > 3. The zero mode solution corresponds to the Ma- 
jorana edge states with Majorana fermions localized at 
the ends of ID wires. To see it clearly, we rewrite the 
representation of the quasi-particle operators in terms of 
Majorana operators: 



vi = 



1 L 



(6) 



where <p nji = (u n ,i-\-v n ,i), r/> n ,i = (u n ,i-v n ,i), if = q+Cj 
and jf = (cj — c\)/i. Here jf and are Majorana 
fermion operators satisfying the relations (7")^ = 7" 
and anticommutation relations {7^,7^ } = 25ij5 a p with 
a and (3 taking A or 5. 

In Fig. 5, we display the distributions of (pi (ipi) for the 
lowest excitation solution of Ai. When V < V c , (pi (ipi) 
is located at the left (right) end and decays very quickly 
away from the left (right) edge. Since there is no overlap 
for the amplitudes of and , the zero mode fermion 
splits into two spatially separated Majorana fermions. In 
the side oiV < V c , the edge mode decays more quickly 
as it deviates farther from the transition point V c . As a 
comparison, distributions of (pi and ipi for the lowest ex- 
citation modes in the regime of V > V c is also given. The 



FIG. 6: Z2 topological invariant vs V for systems with a 
(VE-l)/2, L = 500 and various A. 



solutions of (pi and ipi for V = 3.5 overlap together and lo- 
cate inside of the open chain as a result of Anderson local- 
ization. Consequently, the corresponding quasiparticle is 
a localized fermion which can not be split into two inde- 
pendent Majorana operators. Therefore, the transition 
from TSCs to Anderson localizations can be also judged 
by the presence or absence of Majorana edge fermions in 
different parameter regimes of the system with OBC. 

Z 2 topological invariant.- The existence of Majorana 
edge states is attributed to the nontrivial topological na- 
ture of the bulk superconductor, which can be charac- 
terized by a Z 2 topological invariant [2j . In terms of the 
Majorana fermion operators, the Hamiltonian (1) can be 
represented as the Hamiltonian form 



H 



2L 

E 



A 



A* 



-At. 



(7) 



where L is the number of lattice sites, A is a skew- 
symmetric matrix, 7/ is defined as 72^-1 = jfrfzj = 7 

The nonzero matrix elements 



B 



and {7;,7 m } = 2<5 ; 
are given by A 2j -i t2 j = -A 

A2j-l,2j+2 = —A2j+2,2j-l = A 

-A 2 j+i,2j = 1 + A for j = 1, • • 



_i = Fcos(27rja), 
- 1 and j4.2j,2j+i = 
L with the bound- 



ary condition of L + 1 = 1. For a skew-symmetric matrix 
A, the Pfaffian is defined as follows 



Pf(A) = 



2 L Ll ^ 



sgn(r)A T ( 1)ir ( 2 ) • • • A r (2i-i), T (2i)>(8) 



where S2L is the set of permutations on 2L elements and 
sgn(r) is the sign of permutation. With Pfaffian of a 
system, the Z 2 topological invariant is defined as M = 
sgn(Pf(A)). As shown in Fig. 6, the Z 2 topologically non- 
trivial phase is characterized M — — 1, whereas the Z 2 
topologically trivial phase corresponds to M — 1 . For the 
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system with V < 2 + 2 A, the Z 2 number M — — 1 and 
the system is in the topologically non-trivial phase, while 
for the system with V > 2 + 2A, the Z 2 number M = 1 
and the system is in the topologically trivial phase. As 
the strength of V increases, a topological phase transition 
happens. 

Summary. - In summary, we study the effect of disorder 
produced by the incommensurate potential in ID p-wave 
superconductors which support a topological SC phase 
with Majorana edge states. Increasing the strength of 
disorder destroys the topological SC phase and drives 
the system into a Anderson localized state. The phase 
transition driven by the disorder is identified by analyzing 
the change of gap, the long-range order of the correlation 
function of nonlocal operators and the IPR which char- 
acterizes the spacial localization of wavef unctions. A Z 2 
topological invariant is also used to identify the transi- 
tion from the topological SC phase, which has emergent 
Majorana edge states for the system with OBC, to the 
topologically trivial localized state. 
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